Triplet superconductivity induced by screened phonon 
interactions in ferromagnetic compounds 
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We propose that screened pairing interactions mediated by phonons could give rise to a spin 
triplet superconductivity in ferromagnetic compounds such as UGe2. It is pointed out that the 
pairing interactions include anisotropic components such as those of p, d, /-waves in addition to 
dominant s-wave component due to the momentum dependence. Since the ferromagnetic long- 
range order coexists, there is a large splitting of the Fermi surfaces of up and down spin electrons, 
which suppresses singlet pairing. Therefore, triplet pairing occurs at last due to the sub-dominant 
anisotropic interactions, even in the absence of magnetic contribution to the pairing interactions. 
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Recently, a coexistence of superconductivity and ferro- 
magnetic long-range order was observed in UGe2 under 
pressure In the phase diagram on the pressure and 
temperature plane the superconductivity appears inside 
the area of the ferromagnetic phase. For this proximity 
of the superconductivity and the ferromagnetism, spin 
triplet superconductivity is a possible candidate in this 
compound. 

In addition, singlet pairing is considered to be unfaval- 
able in this compound. The ferromagnetism and super- 
conductivity occur in the same electron band [jl],|| or at 
least in very close electron bands from the crystal struc- 
ture. Thus, there is a large Fermi surface splitting, which 
suppresses antiparallel spin pairing, in the electron band 
responsible to the superconductivity. 

As a mechanism of triplet pairing, magnetically medi- 
ated superconductivity has been considered Q. However, 
there are some behaviors which are not easily explained 
only by this mechanism. For example, the magnetic fluc- 
tuations and their contribution to the pairing interactions 
increase near the transition points. Such behavior is re- 
produced in the calculation of Fay and Appel || . In their 
calculation, it decreases in a narrow region of the width 
of 1% of the exchange interaction parameter near the sec- 
ond order transition point, but this decrease would not 
occur in UGe2 , since the magnetic transition is of first or- 
der in UGe2- Therefore, the superconducting transition 
temperature seems to increase as the magnetic boundary 
is approached in this mechanism. However, in the obser- 
vation, T c decreases near the phase boundary (~ 1.6GPa) 
within the width of ~ 0.3GPa of the pressure, which is 
not very narrow. 

Decreases of the superconducting transition tempera- 
tures as the magnetic phase approached were also ob- 
served in high-T c superconductors. However, in those 
compounds, the magnetic phase is antiferromagnetism. 
Thus, we have a physical explanation of the decrease of 
the superconducting transition temperature based on the 
reduction of the density of states near the Fermi sur- 
face (pseudogap) due to the antiferromagnetic fluctua- 
tions pi. On the other hand, the ferromagnetic fluc- 



tuations do not induce any pseudogap near the Fermi 
surface, since it is split to two surfaces. 

On the other hand, if the decrease of the supercon- 
ducting transition temperature near the magnetic phase 
boundary at the high pressure is correlated to the de- 
crease of the ferromagnetic transition temperature, the 
absence of the superconductivity at pressures p < lGPa 
does not seem to be explained without any extra mecha- 
nisms. For example, pair breaking effect due to Lorentz 
force by internal magnetic field created by ferromagnetic 
moments might suppress the superconducting transition 
temperature at low pressure. It might be also possible 
that the degree of nesting of the Fermi surface g] changes 
by pressure and gives rise to the pressure dependences 
of T c and T x defined in JjJ. Therefore, the magnetic 
mechanisms of the superconductivity is not conclusive 
at present. 

In this paper, we propose a phonon mechanism of 
triplet superconductivity in UGe2 system. In the con- 
ventional mechanisms, the phonon mediated interactions 
have often been considered not to induce anisotropic su- 
perconductivity. However, we shall point out in this pa- 
per that the momentum dependence of the phonon me- 
diated interactions could give rise to a triplet supercon- 
ductivity in ferromagnetic superconductors. In a spin 
polarized state, the spin and charge degrees of freedom 
are locked together, and we concentrate ourselves on clar- 
ifying properties of the phonon mediated interaction. 

It is easily verified that the pairing interactions con- 
tain both singlet and triplet components due to the mo- 
mentum dependence. For example, when the screening 
length becomes longer for weak screening, the pairing in- 
teractions mediated by phonons have shaper peak near 
q ~ in momentum space. Then, anisotropic compo- 
nents have large magnitudes. In the limit where the 
pairing interactions are proportional to the 5-function, 
all components have the same magnitude. Foulkes and 
Gyorffy theoretically examined p-wave pairing in metals 
due to electron-phonon interaction, when the short range 
Coulomb interaction suppresses the s-wave pairing || . 

Anisotropic superconductivity induced by phonon in- 
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teraction has been considered by Abrikosov as a mech- 
anism of high-T c superconductivity In this mecha- 
nism, it was shown that the gap function can vary in 
sign in the presence of on-site Coulomb repulsion. The 
momentum dependence of the pairing interactions was 
examined, where screening effect was taken into account. 
Bouvier and Bok calculated the superconducting gap 
and obtained anisotropic momentum dependence of the 
gap function Ml in the same model. Recently, Friedel 
and Kohmoto gy, and Chang, Friedel, and Kohmoto ||] 
have shown that <i-wave superconductivity is induced by 
screened phonon interactions with an assist of a contri- 
bution from the antifcrromagnetic fluctuations. In the 
studies so far, it was shown that a singlet pairing com- 
ponent is dominant in the screened phonon interactions. 

Usually, since the largest component has even parity 
in the momentum space, it contributes to singlet pair- 
ing. Therefore, a sub-dominant triplet pairing interaction 
cannot induce superconducting transition. However, in 
ferromagnetic superconductors such as UGe2, the dom- 
inant singlet pairing interaction would not be able to 
overcome the Pauli pair breaking effect. Therefore, the 
sub-dominant interaction could give rise to the triplet 
superconductivity. 

In the following, we illustrate that the phonon medi- 
ated pairing interactions contain both singlet and triplet 
pairing interactions by taking into account the screening 
effect. We calculate a coupling constant of sub-dominant 
triplet pairing interactions, which becomes dominant af- 
ter singlet pairing is suppressed. It is shown that the 
screening effect gives rise to an additional lattice con- 
stant dependence of the pairing interactions through the 
screening length scaled by the inverse of the Fermi mo- 
mentum. 

We examine a model with a screened phonon interac- 
tion defined by 



V(k,k') = 



~9— 



(1) 



where q = k — k', q = |q|, and is the screening 

length. This interaction has been examined by many 
authors |6|-||] for the high-T c superconductivity. As ex- 
plained in the text books, near the Fermi surface where 
^ the interaction is attractive due to over- 

screening. We put (£k — £k') 2 ~ for simplicity, since 
k^T -C wo where lod denotes the Debye frequency. 
Hence, we have 



V(k,k') 



|k-k'| 2 



(2) 



Let us first consider the Thomas-Fermi screening for 
a while for a qualitative argument. We will improve our 
theory with a more detailed momentum dependence of 
the dielectric funcion in the random phase approximation 
(RPA) later. In the Thomas-Fermi approximation, we 
put 



47re 2 p t (M), 



(3) 



which is valid for long wave lengths such as q <C kp. 
Here, pt(fJ>) is the total density of states of the conduction 
electrons per unit volume at the chemical potential /i. 

We assume a spherically symmetric Fermi surface for 
simplicity. Assuming top, <C ep, we can put |k| w kp. 
Thus, we have 



v(k,k') 



g(<*-i) 

a — cos 9 



(4) 



where 9 denotes the angle between k and k', and a 
l + q s 2 /2k F . 

The pairing interaction eq. (4) is expanded as 

oo 

V(k,k') = ^^P,(cos#) 

1=0 
oo / 



(5) 



= 9nim(9,fhi m (9',tp'), 



1=0 m=-l 



where (9, <p) and (9',ip') denote the directions of k and 
k', respectively. Here, gi = Vj(2l + 1) and ^yi m (0,(p) — 
Yi m (9,ip) with the spherical harmonic functions Yz m . 
The s-wave coupling constant is calculated as 



go = Vo = 



\g{a - 1) log | — 
2 a 



(6) 
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The coupling constant go is larger than gi, but it is not 
effective in practice, since singlet pairing is suppressed as 
we explained above. Thus, we calculate the next domi- 
nant p-wave component and obtain 



9i = \vi= f V \6) cos 6^- 



g{a-l) 



a . a — 1 

— log 

2 a + 1 



(7) 



In order to improve our approximation, we examine 
the screening effect in an RPA. The pairing interaction 
is modifed as 



V(k, k') 



with 



~ 9 q 2 + q s *u(q/2k F ) w ( q ) 2 - (& _ £ k ,)2 



u(x) 



l r 

2 



f 



2.r 



log 



1+JE 

I -x 



(8) 



(9) 



according to the result of the dielectric function by Lind- 
hard 



«(ff) = 1 



(10) 



The Thomas-Fermi approximation could be recovered in 
the long wave length limit. We apply the same simplifi- 
cation putting (£k — £k') 2 ~ in eq. (8). The coupling 
constant is obtained by numerical calculation. 
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Figure [T] shows the dimensionless coupling constants 
defined by X p = |^i |iV(0), where N(0) is the density of 
states per a spin and unit volume, in the Thomas-Fermi 
approximation and an RPA. It is found that they exhibit 
a peak around q 2 j1k\ ~ 1/2, and the Thomas- Fermi 
approximation agrees well with the RPA for q 2 j2k\ < 
1/2. 
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FIG. 1. Dimensionless coupling constant X p of p-wave pair- 
ing in the unit of gN(0) as a function of q s 2 /2kp. Solid and 
dashed lines show results in the RPA and the Thomas-Fermi 
approximation, respectively. 

The transition temperature of the p-wave superconduc- 
tivity is given by 



T cp = 1.13wDexp[-l/A p ] 



(11) 



in the weak coupling limit. Short range Coulomb repul- 
sion is not very effective in the case of p-wave pairing || . 
Figure || shows the transition temperature of the triplet 
superconductivity based on the formula (11). Because of 
the singular exponential form of eq. (11), the peak of T c 
is sharper than that of X p . 

The screening length q~ 1 is considered as follows. The 
total density of states per unit volume is expressed as 



Pt(e F ) 



m 



Inserting eq. (12) into eq. (3), we have 
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2fcp 7TOH 



(12) 



(13) 



where an denotes the Bohr radius (an = h 2 /e 2 m w 
0.5293A). Since k-p oc 1/a with lattice constant a, this 
expression means that a large lattice constant results in 
a strong screening effect, (i.e., a short screening length 
qs- 1 ). 



8 N(0) = 1. 
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FIG. 2. Superconducting transition temperature T cp of 
p-wave pairing as a function of q s 2 /2kp calculated by eq. (11). 
ojd = 300[K] and gN(0) = 1.6 are assumed as an exam- 
ple. Solid and dashed lines show results in the RPA and the 
Thomas-Fermi approximation, respectively. 

In our spherically symmetric system, the Fermi mo- 
mentum fcp is expressed as k-p = (37r 2 ) 1 / 3 n 1 / 3 /a, where 
n denotes conduction electron (or hole) number per a 
site. Thus, eq. (13) is written as 



q s - I a _ a 

2fcf = n{3n 2 ) 1 / 3 n 1 / 3 m = ai/ 

Here, we have defined a characteristic length ao by 



(14) 



a 



(15) 



which scales the lattice constant. If we put n ~ 1 as 
the order of magnitude, we obtain ao = 2.57A. When 
a — 4A 0,0, one has q 2 j2k\ = a/a — 1.6 as the 
order of magnitude. 

Now, we discuss the pressure dependence of the tran- 
sition temperature of the triplet supercoductivity. The 
pressure dependence of the transition temperature is not 
easily figured out theoretically in the phonon mechanisms 
even in s-wave superconductors. The Debye frequecy wd 
increases with pressure. Then, the prefactor of the weak 
coupling expression of T c increases. On the other hand, 
the electron-phonon coupling constant and the electron 
density of states would decrease. These effects would con- 
tribute to the the pressure dependence in conventional 
superconductors. In UGe2, however, the superconduct- 
ing transition temperature increases more sensitively to 
the pressure near lGPa jt) than those in conventional 
superconductors. 
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FIG. 3. Superconducting transition temperature T cp of 
p-wave pairing as a function of ao/a calculated by eq. (11). 
uj-d = 300[K] and gN(0) — 1.6 are assumed as an exam- 
ple. Solid and dashed lines show results in the RPA and the 
Thomas-Fermi approximation, respectively. 

In the present mechanism, we have an additional con- 
tribution to the pressure dependence through the screen- 
ing effect. At high pressures, the lattice constant a is 
shorten. Then, from eq. (14), the value of q s 2 /2kp is re- 
duced by the pressure. Figure || shows dependence of the 
superconducting transition temperature on the inverse of 
the lattice constant. A large value of 1/a corresponds to 
a high pressure. Since the pressure dependences of wd, 
g, N(0) are not considered, this figure only shows the 
additional contribution though the screening effect. It 
is found that the superconductivity does not occur for 
a > ao and that the pressure would enhance the tran- 
sition temperature of the triplet superconductivity near 
a ~ ao. This result is consistent with the behavior of 
the transition temperature in the experimental phase di- 
agram near lGPa qualitatively, although it is more sen- 
sitive to the pressure in the experiment. 

In addition, inside the region of the ferromagnetic 
phase, Lorentz force created by internal magnetic field 
should suppress pair formation. The superconducting 
transition temperature should be reduced from the value 
that we estimated in the above. Thus, the superconduc- 
tivity must be absent in a wider region than those in 
Fig. H and Fig. [|. 

On the other hand, decreases of the superconducting 
and the ferromagnetic transition temperatures are ob- 
served at high pressures (~ 1.6GPa) Q. For such high 
pressures, crystal structure might become unstable. Sin- 
gle phases might not be able to exist in this region. The 
decreases of the transition temperatures might be due to 
the instability of the crystal. At very high pressures, an 
unharmonicity may play some role in the pressure depen- 
dence of the superconducting transition temperature. 

In addition, near the phase boundary between the fer- 



romagnetic and paramagnetic states, the ferromagnetic 
fluctuations should be strong. The superconductivity 
might be suppressed near the phase boundary by some 
renormalization effect due to the strong magnetic fluctu- 
ations, since the pairing interactions would not increase 
significantly there in phonon mechanisms. In a spin po- 
larized state the spin and charge degrees of freedom are 
locked together. The pressure dependence of T c and T x 
might be explained by taking into account the change in 
the degree of nesting of the quasi-two-dimensional parts 
of the Fermi surface Q . 

In conclusion, we proposed a phonon mechanism of 
triplet superconductivity in ferromagnetic systems. The 
momentum dependence of the pairing interactions me- 
diated by phonons give rise to both singlet and triplet 
components. The singlet interaction is stronger than the 
triplet one, but the Fermi surface splitting suppresses the 
singlet superconductivity in the ferromagnetic state. For 
an appropriate region of the screening length, the triplet 
superconductivity occurs. The superconducting transi- 
tion temperature depends on the lattice constant through 
the screening length scaled by inverse of the Fermi mo- 
mentum. Such dependence might contribute to the pres- 
sure dependence of the superconducting transition tem- 
perature in UGe2- 
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